Elastic Moduli in Nano-Size Samples of Amorphous Solids: System Size Dependence 
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This Letter is motivated by some recent experiments on pan-cake shaped nano-samples of metallic 
glass that indicate a decline in the measured shear modulus upon decreasing the sample radius. 
Similar measurements on crystalline samples of the same dimensions showed a much more modest 
change. In this Letter we offer a theory of this phenomenon; we argue that such results are generically 
expected for any amorphous solid, with the main effect being related to the increased contribution 
of surfaces with respect to bulk when the samples get smaller. We employ exact relations between 
the shear modulus and the eigenvalues of the system's Hessian matrix to explore the role of surface 
modes in affecting the elastic moduli. 



Motivation: This Letter is motivated by some experi- 
mental measurements of the elastic response of nanosam- 
ples of metallic glasses when the radius of cylindrically 
shaped samples was reduced [1]. In this Letter we pro- 
vide a theoretical explanation of this phenomenon. Since 
the experiments are performed at temperatures that are 
much lower than the glass- transit ion temperature, we can 
in the following disregard thermal effect, and study the 
phenomenon in athermal conditions. 

Mathematical background: For the sake of con- 
creteness we study theoretically the shear modulus, 
stressing the difference in its exact calculation between a 
perfect crystalline sample and an amorphous solid sam- 
ple of the same physical dimension. In both cases the 
calculation of the shear modulus starts with the poten- 
tial energy which for a shear-strained solid can be written 
as U({ri (7a/?)}, 7a/?) where {r^}^ are the positions of 
the TV particles and 7^ is the applied strain. (When 
possible we treat j a p as a scalar 7; all the equations can 
be written in tensorial form if required). In both cases 
we can consider the deformation under an infinitesimal 
shear strain via the parameterized transformation on the 
particles coordinates J (7) = I + jh where h determines 
the imposed deformation. Here comes the important dif- 
ference between crystalline and amorphous solids [2] . For 
the former 
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in a perfect crystalline solid (1) 



FIG. 2: (Color online) The dependence of the shear moduli 
\i xz and \i xy on the ratio of the edge size over the height for 
systems varying by a factor of 3.5 in this ratio. One observes 
a change of about 50% and about 15% in fi xy . 



since the particles remain in mechanical equilibrium also 
after the deformation, with all the forces fi vanishing on 
each and every particle. In an amorphous solid, on the 
other hand, even an infinitesimal deformation gives rise 
to non-zero forces on the particles, resulting in non-affine 
displacements Ui that are necessary to restore mechanical 
equilibrium. Thus for an amorphous solid 




FIG. 1: (Color online) A cuboid shape. The top and bot- 
tom surfaces (green) were clamped to measure fi xz , and the 
opposite side wall (red) for \i xy . 



Vi — >• J • Vi + Ui , in an amorphous solid. (2) 

This crucial difference translates to a different calculation 
of the mechanical moduli; in a perfect crystal we can 
write d/dj = d/dj. In an amorphous solid on the other 
hand 
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Applying this to the definition of the shear modulus, 
H = (l/V)d 2 U/dj 2 (with V being the volume) we find a 
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FIG. 3: (Color online) Same as Fig. 2, but with periodic 
boundary condition. 
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FIG. 4: (Color online) The seven lowest eigenvalues of the 
Hessian matrix multiplied by the volume of the sample, in- 
cluding the zero eigenvalue which relates to the Goldstone 
mode. Note the strong decline as Lx/Lz approaches 1. 



different answer, 
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, amorphous solid .(5) 



The shear modulus of the crystalline solid, which is the 
same as the first term in the shear modulus of the amor- 
phous solid is known as the Born approximation. This 
approximation is corrected by the second term which is 
due to the non-affine response of the amorphous solid. 

Eq. (5) is brought to final form using the mechanical 
equilibrium condition dfi/dry = 0, which, using Eq. (3) 
becomes dfi/dj+ (dfi/drj) ■ {duj/dry) =0 [3]. Identify- 
ing the force on the ith particle with fi = —dU/dri, we 
invert the last relation in favor of duj/dj to write 
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Using Eq. (6) in Eq. (5) we get the final result 
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Noticing that the Hessian matrix H is real and sym- 
metric, the difference between the shear modulus of a 
crystalline and amorphous solid is negative definite, nec- 
essarily reducing the shear modulus in the case of the 
amorphous solid compared to the crystalline counterpart 
(with the same inter-particle potential). To understand 
the experimental observations of the mechanical soften- 
ing of smaller and smaller samples we need to understand 



why the correction term increases in absolute magnitude 
compared to the Born term which is system-size indepen- 
dent to a very good approximation. We start by comput- 
ing Eq. (8) using numerical simulations. 

Numerical Simulations: Glassy amorphous sam- 
ples were achieved by using a binary mixture of point 
particles interacting via modified Lennard- Jones poten- 
tial with three different characteristic interaction lengths 
a ss = 1 and <ju = 1.4 and <j s £ = 1.18. Details of the po- 
tentials can be found for example in Ref. [5, 7]. Cuboid 
samples of fixed height (w 15 particles) and square cross- 
section with varying edges (from 50 to 15 particles) were 
prepared by quenching from the melt using a gradient 
energy method to cool the system to T = with zero 
pressure. The boundary conditions were free on all the 
edges except the two edges that were clamped to pro- 
duce infinitesimal shear strain. Clamping the upper and 
lower edges we obtain a strain j xz to measure fi xz while 
clamping two opposite side walls resulted in measuring 
ji xy (cf. fig. 1). For any given system the Born term was 
computed directly from the partial derivatives [4, 8]. The 
correction term was obtained by computing the Hessian 
matrix and the vector H, again directly from their def- 
initions, inverting the Hessian we compute S • H~ x • S 
exactly. Results for both shear moduli as a function of 
the ratio of width to height are shown in Fig. 2. 

We observe a much larger change in \i xz (about 50%) 
than in \i xy (about 15%). We will argue below that this 
difference stems from the contribution of free surfaces. 
In order to emphasize the role of the free surfaces, we 
performed the same measurements of the shear modulus 
in an infinite system made of finite cells which repeat by 
applying periodic boundary conditions. In this case the 
shear modulus exhibits the value of the bulk modulus for 
both directions and for various cell size, cf. Fig. 3. 
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FIG. 5: (Color online) Upper panel: an xz projection by av- 
eraging over the y position of the magnitude of the elements 
of the first nontrivial eigenfunction of the Hessian for the case 
L x /L z — 1.86. Lower panel: the same but for an xy projec- 
tion, averaging over z. The color code is given on the right. 
Note the strong concentration on the surfaces 



Theoretical explanation: Having observed a ten- 
dency towards softening that is very comparable in mag- 
nitude to the experimental one, we are now in a position 
to provide a theoretical explanation. Qualitatively speak- 
ing we expect that particles close to the free surfaces 
should have softer fluctuations in their positions, being 
less constrained than particles in the bulk. By lowering 
the volume of the samples we make the contribution of 
particles at and near the surface more dominant, result- 
ing in an overall softening of the mechanical response. 
Being equipped with the theory presented above, we can 
prove this qualitative intuition by examining the eigen- 
values and eigenfunctions of the Hessian matrix H and 
how they depend on the volume. The negative definite 
term in the amorphous shear modulus can written in a 
way that makes the contribution of the eigenvalues of 
the Hessian explicit by expanding H in the eigenfunc- 
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FIG. 6: (Color online) The same as in the previous figure, but 
for the second nontrivial eigenfunction of the Hessian. Again 
we see that the lowest eigenvalues belong to eigenfunctions 
that are strongly concentrated on the surfaces. 

tions \l> (/c) of H: 

E = £)a fc *< fc >; a fc = E.*(*>. (9) 

k 

With this in mind we can write the correction term to 
the Born approximation in the form 

~v aH — (10) 

where is the eigenvalue of the Hessian matrix associ- 
ated with the eigenfunction ty( k \ 

In Ref. [6] it was shown that \ak\ 2 is roughly system 
size independent. This is not the case for the eigenval- 
ues. In Fig. 4 we show the seven lowest eigenvalues of the 
Hessian matrix (including the Goldstone mode [2]) mul- 
tiplied by the volume, and observe the strong decline in 
these products when L x /L z becomes smaller than, say, 
1.6. This decline will increase the negative correction 
term and accordingly will reduce the moduli. It is also 




FIG. 7: (Color online) An example of a crystalline solid in 
three dimensions. The order of the crystal is distorted due 
to free surface effect. Note that the blue particles are on the 
surface and as the color code tends to red we observe bulk 
crystalline order. 



obvious why \i xz is much more sensitive to this effect 
than fi xy - in the former case, the large top and bottom 
surfaces are clamped and only the side walls are free to 
contribute. Upon increasing the ratio Lx/Lz the con- 
tribution of the free surfaces decreases, and the shear 
modulus increases until we reach the value of the bulk. 
In the latter case, the large top and bottom surfaces are 
free, contributing to the softening of the system. Clearly, 
we need large free surfaces to observe a sizeable effect; 
when we reduce the volume the ratio of surface to vol- 
ume increases, lowering the overall shear modulus. We 
can directly prove that the effect is strongly connected 
to the free surfaces by examining the eigenfunctions as- 
sociate with the lowest lying eigenvalues. In Figs. 5 and 
6 we present the first and second eigenfunction by show- 
ing the magnitude of elements of the eigenfunctions in 
real space. 

We see that both eigenfunctions are strongly concen- 
trated on the surfaces. This is a direct evidence that 
the lowest eigenvalues of the Hessian, which are respon- 
sible for the softest mechanical response, indeed belong 
to modes that live on the surfaces of our samples as pro- 
posed above. 

At this point it is interesting to compare the simula- 
tion results to similar simulations on crystalline samples 
of the same range of sizes and the same geometry. We 
attempted to produce perfect crystalline samples by ar- 
ranging either 'small' or 'large' particles on a lattice at 
zero temperature. We find that when we minimize the 
energy, surface tension effects destroy the crystalline or- 
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FIG. 8: (Color online) The dependence of the shear moduli 
on the ratio of the edge size over the height, for an amorphous 
solid (stars), and a crystalline solid (circles) contains either 
small or large particles, normalized to the value of the shear 
modulus of a cube. 
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FIG. 9: (Color online)The same as in Fig. 4 but for a crys- 
talline solid. Here the decline is weaker than the case of amor- 
phous solid. 



der near the interfaces, leaving us with bulk crystals with 
amorphous interfaces, as seen for example in Fig. 7. 

We thus expect to still see an effect of softening when 
the surface to bulk ratio is reduced, but the effect should 
be smaller compared to the completely amorphous sam- 
ples. For the sake of comparison we show normalized val- 
ues of n xz (normalized to the lowest value) for the three 
types of samples, amorphous and crystalline with small 
or large particles. In Fig. 8 we see that the crystalline 
samples show very similar curves for ji XZl and the effect is 
considerably smaller than for the amorphous solid. The 
relative decrease in the size of the effect is also born out 
by the values of the lowest eigenvalues of the Hessian ma- 
trix (multiplied by the volume) as seen in Fig. 9. The 
eigenvalues times the volume still have a reduction for 
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the largest values of the surface to volume ratio, but the 
degree of reduction is considerably smaller, cf. Fig. 4. 
We expect that for larger samples (as in the experiment) 
where the surface to bulk ratio is smaller than in our 
simulations, the effect will be even smaller than in our 
present simulations. 

In summary, we have shown that the main reason that 
is responsible for the softening of the mechanical response 
of nano-samples as a function of their size comes from the 
softer normal modes of the hessian matrix that are con- 
centrated near the boundaries. These reflect the physics 
of increased freedom of particles near the surfaces com- 
pared to the highly constrained particles in the bulk. 
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